function [act logact] = gmmactivBayes(mix, data)


N = size(data.Y, 1);
K = mix.ncentres;

m = mix.varposterior.m;
b = mix.varposterior.b;
W = mix.varposterior.W;
v = mix.varposterior.v;

act = zeros(N, K);
logact = zeros(N, K);
for k = 1:K
    [act(:, k) logact(:, k)] = gaussBayes(m(k, :), b(k), W(:, :, k), v(k), data.Y);
end
